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Abstract 

We look for the necessary conditions allowing the Universe isotropi- 
sation in presence of a minimally coupled and massive scalar field with a 
perfect fluid. We conclude that it arises only when the Universe is scalar 
field dominated, leading to flat spacelike sections and accelerated expan- 
sion, and examine the case of a SUGRA theory. 
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1 Introduction 

Today the introduction of scalar fields in cosmology obeys major reasons taking 
roots in the extension of the standard particle physics model such as super- 
symmetry which requires additional degrees of freedom represented by these 
fields. In a general way, most of the theories predicting extra dimensions at 
high energy could generate scalar fields^ via compactification processes. Su- 
pergravity theory (SUGRA) El related to supersymmetry concept or Higgs 
mechanism which allows us to explain the mass of particles also imply some 
scalar fields. From an observational point of view they could be responsible for 
dark matterj^I E] as well as dark energy[51 ^ ^1 [7] although other explana- 
tions exist. They could also solve the so-called cosmological constant problem: 
most of these scalar fields are massive and thus able to mimic a variable cosmo- 
logical constant. 

Let us speak about the geometrical context of this paper. There exist nine 
anisotropic cosmological models classified by Bianchi in 1897. We will be inter- 
ested in the curved Bianchi class A models since we have studied the spatially 
flat Bianchi type / model in ^I] and there is no adapted ADM Hamiltonian 
formulation for the Bianchi class B models. Among the Bianchi class A models, 
the Bianchi type IX one contains the solutions of the positively curved isotropic 
FLRW model and the Bianchi type // one characterizes the strong anisotropic 
phases m. Unless we assume a Universe born isotropic and homogeneous, as 
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instance thanks to a quantum principle selecting this type of particular model 
among all the possible ones, it is legitimate to ask why our Universe is so sym- 
metric. It seems more natural to suppose that it was initially less symmetric and 
that it asymptotically evolves to an FLRW model. This is one of the reasons 
why the study of anisotropic models is so important. It allows us to explore 
the mechanisms responsible for the isotropisation of our Universe and to put 
some constraints that may be compared to observations on its final isotropic 
state. Moreover, from the initial state point of view, the oscillatory approach of 
the singularity by the Bianchi type IX model is generally considered as more 
generic than the one of the FLRW models and could be shared by the most 
general inhomogeneous models as conjectured by Belinskij, Khalatnikov and 
Lifchitz pi[Tl| . 

Our goal is to find some scalar field properties allowing the Universe to reach 
isotropy and then the dynamical behaviours of the metric and potential. In |15j , 
we have shown that isotropisation of curved class A Bianchi models in presence 
of a massive scalar field but without a perfect fluid always leads to a late times 
acceleration which is not necessary the case when there is no curvature|16j. In 
[TT] , we have seen that in presence of a perfect fluid, the isotropisation of the 
flat Bianchi type / model leads to a decelerated expansion if asymptotically the 
difference — between the pressure and the density of the scalar field is pro- 
portional to the density p of the perfect fluid. What happens when we consider 
both curvature and perfect fluid? Here, we will try to answer this question. 
To this end, we will use the ADM Hamiltonian formalism If to get a first order 
equations system that we will study by help of dynamical systems analvsis|18|. 
Most of times, dynamical analysis of the field equations in cosmology rest on 
the orthonormal frame formalism and Hubble-normalized variables as shown in 
Wainwright and Ellis book^H]- It allows us to study a large number of cos- 
mological models in various situations, even the most complex one such as the 
inhomogeneous cosmologies'!^ or the presence of magnetic fields [SU], finding 
and classifying all the equilibrium points of these systems. Some scalar-tensor 
theories have also been studied in this way but, to our knowledge, their forms 
were always completely specified, i.e. they did not contain any unspecified func- 
tion of the scalar field. Here, we want to consider a class of scalar-tensor theories 
containing two unspecified functions of the scalar field and just look for the sta- 
ble isotropic state the Universe can reach. Hence, we aim to study a larger 
class of scalar-tensor theories than usually and it is one of the reasons why we 
have not used the powerful orthonormal frame formalism but rather the more 
traditional Hamiltonian ADM formalism which have proved to be useful in such 
a case[?T1. 

The plan of this work is as follows: in the second part we establish the Hamil- 
tonian field equations and, after having remembered the results we obtained 
without a perfect fluid, we study the isotropisation process when it is present. 
We discuss the physical meaning of our results in the last section. 

2 Field equations and dynamical analysis 

In the first subsection, we derive the Hamiltonian field equations and in the 
second one, we use dynamical systems analysis to study the stable isotropic 
states. 
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2.1 Field equations 



We will use the following metric, reflecting the 3+1 decomposition of spacetime: 

ds'^ = -{N^ - NiN')dn^ + 2NidVLuj' + Rlg.juo^uj'^ (1) 

The uji are the 1-forms generating the Bianchi homogeneous spaces, N and Ni 
are the lapse and shift functions and gij are the metric functions parameterised 
by Misner|22| as: 

g-2n+/3+ + \/3/3_ 
g-2n-2/3+ 

The (3± functions describe the anisotropy whereas is the metric isotropic part. 
The action of the minimally coupled and massive scalar field theory with a non 
tilted perfect fluid writes: 

S^{16tt)-^J [R-{3/2 + uj{(l)))(t)-^'(l)^f,(t)-^-U{(j)) + 16TTc^L^]y/^d*x (2) 

U is the potential of the scalar field <j) whose coupling with the metric is described 
by the Brans-Dicke coupling function the Lagrangian of the non tilted 

perfect fluid whose equation of state is p = (7— l)p with 76 [1,2]. It describes a 
dust fluid when 7 = 1 and a radiative fluid when 7 — 4/3. The other important 
value is 7 = and corresponds to a cosmological constant which has been 
discussed in JS] . Technical details allowing to get the ADM Hamiltonian from 
the action (0) have been given in [JUI^E]. Hence we write directly: 

H^=p\+ pi + 12^^ + 2AT:^Rle-^''U + (Je^^^-^)^^ + V{n, [3+, (3.) (3) 
3 -\- 2u) 

with p± and p^, respectively the conjugate momenta of the (3± variables and the 
scalar field. V^(ri, /?+, /3_) is the curvature potential characterising each curved 
Bianchi class A model and given in tabled 5 is a positive constant proportional 



Type 


Expression of V{n, (3+, /3^) 


II 


247r2i?4e-4f^+4/3++v'3/3- 










VIo 


247r2i?4e-4n+4/3+ ^^^^^ 4 V3/3_ 


+ 1) 








VIIo 


247r2i?4e-40+4/3+ (^^gj^ 


-1) 








VIII 


247r2i?4e-4n [g4/3+ ^^^^^ 


-1)- 


hl/2e-8'3+ 


+ 2e-2'5+ cosh 


-2V3/3-] 


IX 


247r2i?4e-4n [g4/3+ ^^^^^ 


-1)- 


hl/2e-8/3+ 


- 2e"2'3+ cosh 


-2\/3/?_] 



511 = 
522 = 

533 = 



Table 1: Curvature potentials for Bianchi type //, VI^-, VIIq, VIII and IX 
models 

to (7 — l)po- Using (O, the Hamiltonian equations are: 



scalar field transformation sometimes allows to reduce the two unspecified functions uj 
and [7 to a single function. However, tfie transformation is not always analytically possible 
and it is why it is more general to consider the two functions. 
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9^ " (3 + 2uj)H 



(5) 



dH dV 

In this paper, we will choose Ni = 0, i.e. a diagonal metric, and we derive [T7| 
that 

H 

Now, we have to rewrite these equations with some variables, bounded in the 
neighbourhood of the isotropy. In [TS] we had used the following variables 
common to all the curved Bianchi class A models: 

x±=p±H-' (9) 

y = T:RlVUe-^^^H-^ (10) 

z=p^<j>{3 + 2ujy^/^H~'^ (11) 

and (f) the scalar field. These variables are real as long as C/ > and 3 + 2cj > 
which is necessary to respect the weak equivalence principle. Each of them has 
a physical interpretation: 

• are proportional to the shear parameters S± defined in jl8|. 



y"^ is proportional to (p^ — p^) / {dil / dt)^ , (dil./dt)'^ being the Hubble vari- 
able when the Universe is isotropic, and p^ the density and pressure of 
the scalar field. 

• is proportional to {p^+^) / {d^l/ dt^ , (dQ/dt)^. 

• We deduce from these two last points that the density parameter ft^ for 
the scalar field is a linear combination of y'^ and or, when the scalar 
field is quintessent, that these two variables are proportional to fl^. 

We had also defined some "w" variables characterising the curvature of each 
Bianchi model and which are shown in the table|21 They are related to the three 
Ni variables describing the curvature in the paper of Horwood and Wainwright ^ 
or in the book edited by Wainwright and Ellis JHI and defined by using a sym- 
metry group structure. In this last book, the curvature of the Bianchi type // 
model, VIq and VIIq models, VIII and IX models are respectively described 
by A^i , {N2tN^) and {Ni, N2, N3) variables. Here, in a similar way, we could 
redefine three variables Wi, i = 1,2,3 such that for the Bianchi type // model, 
VIo and VIIq models, VIII and IX models, the curvature be described by 
{wi — w), {wi — w+,W2 — W-) and {wi — WpW^,W2 — Wp/w^,w^ = w„i), thus 
recovering the same unified picture as in (IBj . 

In this paper, we will also consider an additional variable called k and related 
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Bianchi models 


Associated variables 


// 




VIo and VIIq 


W± = ^i?2g-2n+2(/3+±V3/3-)^-l 


VIII and IX 





Table 2: w variables characterising he curvature of each Bianchi model, 
to the presence of a perfect fluid. It is defined by 

fc2 = Sy^V-'U-^ (12) 

k is proportional to the density parameter of the perfect fluid, one of the main 
parameters in cosmology. It can be shown by checking that cx V^^ /{^)^ . 
k is not independent from the other variables and when no perfect fluid is 
considered, fc = strictly. 

For each Bianchi model, we have rewritten the Hamiltonian constraint and the 
field equations with these variables in the second appendix. 

2.2 Isotropisation 

In the first subsection, we define the different ways to reach a stable isotropic 
state. In the second one, we recall our results obtained without the perfect fluid. 
In the third one, we discuss about their stability. In the fourth one, we extend 
them by considering the presence of a perfect fluid. 

2.2.1 Different kinds of isotropisation 

In when no perfect fluid is present, we had defined the isotropy as the 
convergence of the metric functions to a common form such as the Hubble pa- 
rameter is the same in any directions. It implied d(3±/dt and P± const 
and thus that it should arise when p±e^^ — > 0. This definition is unchanged in 
presence of a perfect fluid. 

Different kind of isotropisation may exist, leading to a forever expanding model, 
a singularity or a static Universe. We had shown that when there is no perfect 
fluid, isotropy only occurs when — )■ —oo and x 0, i.e. for a forever expand- 
ing Universe. Looking at the field equations, we find three ways to reach an 
isotropic stable state that we have classified in three classes: 

1. Class 1: all the variables but not necessarily the scalar field reach equilib- 
rium with y ^ 0. 

2. Class 2: all the variables but not necessarily the scalar field reach equilib- 
rium with y ~ 0. 

3. Class 3: all the variables do not reach equilibrium but x± which, as the 
w functions, have to vanish. 
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For the class 2, generally nothing can be deduced about the asymptotic be- 
haviours of the metric functions and potential. It has been numerically observed 
in a paper in preparation where a non minimally coupling between a scalar field 
and a perfect fluid is considered. For the class 3, y and z do not necessarily 
need to reach equilibrium when the Universe isotropises. They just have to 
be bounded when ~f — oo, implying that they oscillate. Hence, the signs of 
their derivatives, which do not asymptotically vanish^, change continuously^. 
We have numerically observed the class 3 isotropisation in presence of several 
scalar fieldsjSS) and it seems to be associated to an oscillating behaviour of £. 
In this paper, we will study the class 1 isotropisation. In the two next subsec- 
tions, we briefly recall the results we obtained in jl5| without a perfect fluid and 
then discuss the assumptions we have made related to their stability. 

2.2.2 Without the perfect fluid 

The results we obtained in JHI are the following. We define the function i of 
the scalar field: 

£ = 0C/^[/-i(3-t-2w)~i/2 

The equilibrium points corresponding to class 1 isotropisation are given by 
{x±,y,z) = (0, ±^3 — £2(6-\/2)~-^, ^/6), the w variables related to the curva- 
ture (see table ^ being zero. Our conclusion about isotropisation, valid what- 
ever the curved Bianchi class A models when no perfect fluid is present, was 
that it occurs for a forever expanding Universe {Q, —> ~oo) and it requires £^ 
to tend to a constant £o smaller than 1. Then the metric functions tend to 
t^o if £q ^ or to a De Sitter model otherwise. The Universe is thus asymp- 
totically accelerated and flat. The scalar field asymptotical behaviour may be 
determined by the asymptotical solution of the first degree differential equation 
= 202[/^(3 + 2u;)-i[/-i. 

For Bianchi type //, VIq and VIIq model, isotropisation will occur at late times 
if the Hamiltonian H is initially positive and at early times otherwise. It is eas- 
ily shown by noting that H is a monotonic function of f2 with a constant sign. 
Then, using the relation dt = —Ndfl, it comes that J7 is a decreasing(increasing) 
function of the proper time t when H is positive (negative). Since the Universe 
only isotropises in Q —^ — cxd, it thus corresponds to late times and forever 
expanding Universe. For the Bianchi types VIII and IX models, it is not pos- 
sible to show that is a monotonic function and thus, the isotropisation time 
is undetermined. 

2.2.3 Stability of our results 

The above results or the ones of the present paper are the determination of the 
isotropic equilibrium points, some necessary conditions for isotropisation and 
the asymptotical behaviours of some functions in the neighbourhood of these 
points. However the asymptotical behaviours are determined by calculating the 
exact solutions for each equilibrium point and they will be correct only if on 
one hand £ and in the other hand the variables {y,z,w) (and k when we will 

^We are assuming that they do not reach equihbrium! 
^The variables are bounded 

*It is necessary that £ tends to a constant otherwise, z as instance, which tends to 
could not be asymptotically vanishing. 
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consider a perfect fluid) tend suflicicntly fast to their equilibrium values. Oth- 
erwise, they will be different. Let us explain why. 

The first kind of instability conies from L As instance, when we look for x± 
asymptotical behaviours, we need to calculate exp{J i^dfl). As explains in the 
section VZ.'Z.'ZX we have shown in [TS] that near the isotropic state, P tends to 
a constant £o (vanishing or not). Then, in our calculation, we have assumed 
that asymptotically when Vl — s- —oo, exp{J I'^dft) — > exp{£QQ) but this is true 
only if tends sufficiently fast to its constant equilibrium value. As instance, if 
^ £l + Q^^/^, £^ tends to a constant but exp{J i'^dfl) does not tend asymp- 
totically to an exponential because of the quantity Hence, our results 
will be valid as long as the following assumption holds: 

• When £ tends to a constant £o (vanishing or not) such that £'^ ^ £q + S£'^, 
J{£l + S£^)dn £ln + const. 

If the last limit is not true, the asymptotical behaviours for the metric func- 
tions (and potential) are different from classical power or exponential laws. This 
problem could be overcame since our results allow to calculate 4>{^) and thus 
£{n). Hence, it should be easy to generalise them by keeping the / £^dn term 
instead of considering that it tends to £^fl but then they would not be on a 
closed form. 

The second kind of instability can not be solved so easily. In the same way, 
the asymptotical behaviours we have determined will be true only if the vari- 
ables {y,z,w,k) tend sufficiently fast to their equilibrium values. As instance 
near isotropy we have y — > ±-\/3 ~ £'^{6^/2)^^ and when we integrate the dif- 
ferential equation for x±, we assume that exp{J y^dH.) — > exp{J (3 — £'^)/72d^l). 
But once again, this is not exact if y^ tends to its equilibrium value slower than 
r2~^ and we have to make the same kind of assumption for [y,z,w,k) as for 
£. For partly solve this problem, it would be necessary to consider some small 
perturbations of the exact solutions but until now we have not succeed to get 
any interesting results, even for the flat model. 

To summarize, the results of this paper related to asymptotical behaviours will 
be valid for a class 1 isotropisation if the function £ and the variables (y, z, w, k) 
tend sufficiently fast to their equilibrium values or, more physically, if the Uni- 
verse tends sufficiently fast to its isotropic state. The restriction on £ may be 
easily solved but the ones on (y, z, w, k) require a more careful examination. In 
the following subsection, we consider the isotropisation of a curved Bianchi class 
A model in presence of a perfect fluid, first when k vanishes and then when it 
tends to a non vanishing constant. 

2.2.4 With a perfect fluid 

Depending on the vanishing of k near an isotropic equilibrium state, the results 
summarize in the section [2 . 2 . 21 will or will not be modified. 

fc ^ 

When k ^ Q near isotropy, the results are the same as those found when we 



7 



consider no perfect fluid. In particular, isotropisation always arise for a for- 
ever expanding Universe, i.e. when Q — > — oo. Obviously, we find the same 
equilibrium points and assuming that k tends sufficiently fast to its equilibrium 
value(see section l?.2.3|l . we also recover the same asymptotical behaviours. How- 
ever, the limit fc ^ plays the role of a new constraint. This fact was noted 
in for the flat Bianchi type / model. In this last paper we had shown that 
the interval of £ allowing for isotropy was smaller when we consider a perfect 
fluid such that k —> than without it: in this last case isotropy requires £^ < 3, 
otherwise i'^ < 3/27^. 

Does the limit fc ^ also change the necessary conditions for isotropy when 
we consider some curvature? Near equilibrium, the w variables have to vanish 
and are proportional to e~^^H~^. But e^^^^ diverges and thus, since w ^ 0, H 
have to be larger than e^^^, i.e.: 



Moreover, we have = Ue ^^e ^ and since near isotropy y tends to a 

non vanishing constant whereas e~^^H~'^ tends to vanish, we deduce that 

U » e^^^ » V-^ 

and then from (|12|l that fc ^ 0. Consequently, starting from the fact that w 0, 
we conclude that fc ^ without any modification of the necessary condition 
for isotropisation on the contrary from the Bianchi type / model. Moreover, 
it means that the energy density of the scalar field and its pressure are 
such that U <x — >> pm- the Universe is dynamically dominated by the 
scalar field. Thus, the results obtained in the vacuum (i.e. fc = strictly) are 
not changed when we consider a perfect fluid such as fc — s- 0. 

fc 7^ 

Now, we consider what happens when fc 0. The necessary condition for 
isotropy is still p±e^^ — > and we have to determine if it occurs for a forever 
expanding, contracting or static Universe. 

• If it arises for a diverging f2, it means that at equilibrium, we must have 
a;± — > as explained in the section 12.2. II 

• If it arises for a finite value of f2, then we must have p± 0. 

— Let us assume that in the same time x± -/^ 0. Since p± 0, from 
© we deduce that H have to vanish otherwise a;± — > 0. But then fc^, 
which is proportional to the perfect fluid density parameter, diverges 
and the constraint is not respected because, near isotropy, all the 
variables have to be bounded as shown in |15| . Hence, H can not 
tend to zero and x± must vanish near equilibrium. 

— In the same way, when fl tends to a non vanishing constant, H can 
not diverge because then fc — s- which is not in agreement with the 
assumption of this subsection. 

^This inequality rests on the asymptotical behaviour of k and, as discussed in section|T^21 
it may vary if k does no tend sufficiently fast to its equilibrium value. However, the limit £^ < 3 
have always to be respected since it is required for the existence of the equilibrium points, 
independently on how fast the isotropic state is reached. 
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Consequently, when the isotropy occurs for a finite f2, the Hamiltonian 
have to tend to a bounded and non vanishing quantity and it is thus the 
same for the w variables. 

To summarize, in the neighbourhood of the isotropic state: 

• If f7 diverges, the equilibrium points are such that x± — > 0. 

• If r2 ^ const =/= 0, the equilibrium points are such that x± — > and w 
variables are non vanishing and bounded. 

Whatever the curved Bianchi models, the only equilibrium points correspond- 
ing to these requirements when solving the field equations^ are defined by 

{x±,y, z) — (0, ± ; the w variables related to the curvature(see table 

13) being 0^. The Hamiltonian constraint implies that k"^ = 1 — ^ and conse- 
quently £^ > I7. This inequality is independent from any assumption on how 
far the isotropic state is reached. As the w variables are vanishing, it follows 
that r2 must diverge and not tend to a constant. Consequently, we calculate that 
asymptotically the Hamiltonian, whose form is given in the second appendix as 
a function of the variables, behaves as: 

This is in agreement with the limit const ^ and the definition for 

k which also implies that U oc V~'' . Hence, the scalar field plays the same 
dynamical role as the perfect fluid and we can show that their energy densities 
scales in the same way, preventing any accelerated expansion. We also get that 
the w variables (but W- which tends to a non vanishing constant for the Bianchi 
type VIII and IX variables) all behave as: 

For the considered range of 7, we derive that w ^ only if f2 — > +00. But for 
the x± variables, it comes: 

X[) being an integration constant. It follows that if 7 S [1,2] and Vl +00, x 
diverges. Consequently, the isotropic state can not be reached for this range of 
7. Knowing x± and H, we calculate that: 

Hence, p±e^^, the w and x± variables vanish only if 7 < 2/3 and ft — s- —00. 
Then, we find that — > t~ and, from the definition of y and the property 
U (X V~'^ , we derive that U t~^. This restriction on 7 does not exist for the 
flat Bianchi type / modelj^m^HI and does not fit an ordinary perfect fiuid such 
that 7 e [1,2]. 

®For the Bianchi type VIII and IX models where the equations are far from being simple, 
it is not possible to solve them directly. We proceed by putting x± = and ui_ = 1 in 
the equations for x±, these values being these required for isotropy. We then show that x± 
can reach equilibrium only if Wp and Wm vanish which allow us to determine the equilibrium 
values for the other variables. Since Wp and Wm tend to vanish, the Hamiltonian constraint 
shows that all the variables are bounded. For the other Bianchi models, we can show in the 
same way as in 1151 that all the variables are bounded near the isotropic equilibrium state. 

^There exist some other equilibrium points for which k or £ may be chosen such that x± = 
and the constraint be respected but they correspond to complex values of some variables. 
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3 Conclusion 



In this work, we have determined the necessary but not sufficient conditions 
for class 1 isotropisation of curved Bianchi class A models when a minimally 
and massive scalar field with a perfect fluid are considered. We have assumed 
that U > and 3 + 2aj > such that the weak energy principle is respected. 
Moreover, some of our results related to asymptotical behaviours are valid as 
long as the isotropic state is reached sufficiently fast. 

We can distinguish two cases depending on the vanishing of k, a variable 
proportional to the perfect fluid density parameter Qm- When isotropy occurs 
with fc ^ 0, we have thus ilm 0, U (x — pc/, > Pm and the results are the 
same as in pTSj where no perfect fluid is present: 

Class 1 isotropisation with flm — > 0; 

A necessary condition for isotropisation of curved Bianchi class A models in 
presence of a minimally and massive scalar field such that flm — will be that 
the quantity £ = (l)U^U^^{3 + 2lij)^^^^ tends to a constant £q, whose square is 
smaller than one. For the Bianchi type II , VIq and VIIq models, it arises 
at late (early) times if the Hamiltonian is initially positive(negative). For the 
Bianchi type VIII and IX models, the time of isotropisation is undetermined. 
// £o 7^ 0, the metric functions tend to a power law t o and the potential van- 
ishes as If £o = 0, the Universe tend to a De Sitter model and the potential 
to a constant. The isotropisation process always leads to a fiat and accelerated 
Universe. 

Considering the limit near isotropy of the Hamiltonian equation for (f) rewritten 
with the normalised variables (see appendice), we deduce that the scalar fleld 
asymptotically behaves as the limit of the solution for 

4> = 24>'^U^U-\3 + 2uj)-'^ 

as Q ^ —oo, in the same way as in |15| . This last equation allows us to deduce 
the asymptotical behaviour of £{Q) when we specify lo and U. The second result 
of this work concerns the case for which fc, or equivalently f2m, tends to a non 
vanishing constant implying that [/ cx — cx We have then: 

Class 1 isotropisation with f2„i const ^ 0; 

The isotropisation of curved Bianchi class A models in presence of a minimally 
and massive scalar field such that Q,n ^ const ^ is impossible if the perfect 
fluid is an ordinary one such that 7 G [1,2]. It will only occur if "f < 2/3, which 
generally corresponds to a quintessent fluid equation of state. 

Now, we examine these results with respect to supergravity. In [21 E|, it is 
shown that quintessence theories should be based on supergravity. A scalar 
tensor theory is then derived, defined by w + 3/2 = (fp' and U = A^+'"(/)~'"e'S''^ . 
It is able to solve the coincidence problem and even the fine tuning problem if 
TO > 11. What about class 1 isotropisation? We calculate that: 

2 _ , n(j)^ -m 2 
^ ^ ^/20 ^ 
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and when no matter is present or if k 
behaves as: 



0, the scalar field asymptotically 



4>oe 



2nn 



± 



(po being an integration constant. 

When n > 0, (j) ^ (m/nY^'^ implying that m should be positive, i —> and 
the necessary conditions for isotropisation are thus respected. If it arises, the 
Universe tends to a De-Sitter model. It could thus describe the inflationary 
period, before the domination of the matter. This case is plotted on figureQlfor 
the Bianchi type IX model. 

When n < 0, the scalar field behaves as 6 ±\ jt is defined when 




Fig Ure 1] Isotropisation of SUGRA theory for Bianchi type IX modeh 
tial conditions and potential parameters respectively are {x , x - , z , Wp^ Wm , 
(0.05, 0.83, 0.025, 0.12, 0.02, 0.002, 0.2, 0.14) and (A, m, n) = (2.0, 1.1, 0.15). 



Ini- 



51 — oo if 00 > and then diverges. It follows that £ also diverges and thus a 
class 1 isotropisation is not possible as confirmed by numerical simulations. 
Summarising, if the Universe isotropises, this theory issued from SUGRA leads 
an anisotropic curved Universe to a flat isotropic De Sitter one dominated by 
the scalar field. It could be a good description for an inflationary period. Nu- 
merical simulations have not shown any class 2 or 3 isotropisation. 



In conclusion, we knew that when no perfect fluid is present, the class 1 isotropi- 
sation of an anisotropic curved Universe may lead the Universe to flat spacelike 
sections and accelerated expansion if some necessary conditions are respected. 
The question was: does this acceleration, due to the presence of curvature, al- 
ways exist in presence of a perfect fluid. The answer is "y^s" when the density 
parameter of the perfect fluid asymptotically vanishes. Then, its presence does 
not change the asymptotic isotropic state or the necessary conditions to reach 
it. Contrary to the flat Bianchi type / model for which an isotropic state such 
that flrn — > const 7^ may exist, the perfect fluid and the scalar field playing 



11 



the same dynamical role, the isotropic state in presence of curvature is always 
scalar field dominated but if the perfect fluid is an exotic one. Future research 
should be concerned by a scalar field which violates the energy conditions, i.e. 
such that uj < — 3/2 or J7 < or which is not minimally coupled to a per- 
fect fiuid. This last possibility, which would allow to extend our results to the 
Hyperextended Scalar Tensor theory (i.e. with a varying gravitation function) 
with a perfect fluid, is currently under consideration in a paper in preparation. 
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A Field equations of the curved Bianchi models 
with normalised variables 

Bianchi type // 

The Hamiltonian constraint writes: 

xl+x'i+ 24?/ + I2z^ + I2w^ + fc2 = 1 (13) 

The Hamiltonian equations become: 

x+ = 72y^x+ + 24w^x+ ~ 2\%i? - 3/2(7 - 2)k^xj, (14) 

X- = 72/a;_ + 2Wx^ - 24\/3u;2 - 3/2(7 - 2)k^x- (15) 

y = y{Uz + 72/ - 3 + 24w2) - 3/2(7 - 2)k^y (16) 

z = /(72z - 12^) + 2^w^z - 3/2(7 - 2)k^z (17) 

w = 2w(a;+ + %/3a;_ + 12^2 + 36/ - 1) - 3/2(7 - 2)k^w (18) 

To get an autonomous system, we need a first order equation for </). Rewriting 
it comes: 

This equation is the same for any Bianchi models. The equation for H may be 
rewritten as: 

= -ff(72/-|-24t«2 + ^(7-2)fc2) (20) 

Bianchi VIq and Vll^) models 
The Hamiltonian constraint writes: 

x\ + x^_ + 24/ + 12z2 + 12(u;+ ± w-f + fc^ = 1 (21) 

and the Hamiltonian equations become: 

x+ = 12y'^x+ + 24(a;+ - l)(w_ ± w+f - 3/2(7 - 2)k^x+ (22) 
i_ = 72y2a;_ + 24a;_(w_ ± w+)2 + 24\/3(w2 - w\) - 3/2(7 - 2)fc2x_ (23) 
y = y{Uz + 72y2 - 3 + 24(u;_ ± w^f) - 3/2(7 - 2)k^y (24) 
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i = y'^{72z - + 2Az{w- ± w+f - 8/2(7 - 2)k'^z (25) 
w+ = 2w+ [x+ + \/3a;_ + 12(«;_ ± w+f + 36^^ - l] - 3/2(7 - 2)k'^w+{2Q) 
w- = 2w- [x+ - VSx- + 12{w- ± w+)2 + 36y2 - l] - 3/2(7 - 2)k'^w-{27) 



The equation for H is: 

H=-H 



72y^ + 2A{w+ ± W-f + -(7 - 2)fc^ 



(28) 



Biaiiclii VIII and /A' models 
The Hamiltonian constraint writes: 



C+ + a;2 + 2%2 + I2z2 + I2[wl{l + wt) ± 2w^{wmWpf/'^{l + 10"^)+ 
wl{wl, - 2w^)]{wlwp)-'^ = 1 



and the Hamiltonian equations are: 

x+ = 72y'^x+ + 24{w3(x+ - + wt) ± w-{l + 2x+){w„iWpf/'^{l + wl) 
+wl_ [(2 + x+)wl - 2{x+ - l)wl]}{w^_wp)-^ - 3/2(7 - 2)k^x+ (29) 
X- = 72y2x_ + 24{w3 [wt{x- - VS) + a;_ + V3)] ± W-{wmWp)^^'^[w'i 

{-V3 + 2X-) + (73 + 2x_)] + wlx-{wl, - 2w^)}{wlwp)-^ (30) 
-3/2(7 - 2)fc2a;_ 
y = y{Uz + 722/2 _ 3 + 24[ii;3(l + «;!) ± 2(«;„u;p)3/2w_(l + wl.)+ 

wl{wl - 2wl)]{wlwp)-^} - 3/2(7 - 2)khj (31) 
z = 2/2(72^ - 12£) + 24z[«;3(l + ± 2{w^Wpf'^w^{l + wl)+ 

(u,3^ _ 2wl)]{wlwp)-^ - 3/2(7 - 2)A:2z (32) 
Wp = Wp{-2 + 2x+ + 722/2 + 24[w3(l + wt) ± 2w-{wynWpf''^{l + w'^_) 

(u;^ - 2wl)]{wlwp)-^} - 3/2{j-2)ewp (33) 
Wm = M'm{-2 - 2x+ + 72y2 + 24[u;3(i + wi) ± 2w_(u;„Wp)3/2(l + ^2 ) 

+wUwl - 2wl)]{wlwp)-'} - 3/2{^-2)k^Wm (34) 
w- = 2\/3w-X- (35) 



and 

H = -H[72y' + 24(±2^^^ ± 2wl/'wTw. - 2^1 + <_ + 

t«X + ^) + |(7-« (36) 
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